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We discuss the BFKL equation with a running gauge couphng and identify in its 
solutions the contributions originating from different transverse momentum scales. We show 
that for a running couphng constant the distribution of the gluons making up the BFKL 
Pomeron shifts to smaller transverse momenta so that the dominant part of Pomeron can 
have a nonperturbative origin. It is demonstrated how this soft physics enters into the 
BFKL solution through the boundary condition. We consider two kinematical regimes 
leading to different behaviour of the rapidity and transverse momentum dependence of the 
gluon distribution. In the diffusion approximation to the BFKL kernel with running Us , we 
find a sequence of poles which replaces the cut for fixed Oa- The second regime corresponds 
to the singular part of the kernel, which gives the dominant contribution in the limit of 
very large transverse momenta. Finally, a simple more general picture is obtained for the 
QCD Pomeron in hard processes: it is of soft, nonperturbative nature, but has hard ends 
of DGLAP-type. 



1. INTRODUCTION 



The description of high energy reactions in terms of the exchange of a Pomeron and of secondary 
Reggeons continues to be very successful. Ever since the arrival of QCD, the challenge has of course 
been to provide an elementary, microscopic description of the Regge theory. Nonperturbative aspects have 
so far made it impossible to have a complete QCD description. However, in the extreme limit of very 
high rapidities and large transverse momenta, it was possible to model the Pomeron perturbatively by the 
exchange of a gluon ladder and obtain an analytic expression under certain kinematical conditions. This 
was called the BFKL Pomeron ||l[], which was considered to be of somewhat academic interest. The recent 
HERA experiments, which showed a very rapid increase of the structure function F2(x, Q^) in the limit of 
small X, have renewed the interest in the BFKL Pomeron, since it predicts such a behaviour. This inspired 
the hope that this and other features resulting from the BFKL Pomeron already start to appear at existing 
energies. 

It is commonly assumed that the Pomeron in QCD is dominated by gluons. Due to its quantum 
numbers it has to consist of at least two gluons, which would result in a constant cross section as a function 
of energy The admixture of multigluon states and other, nonperturbative effects then should lead to 
the Pomeron, with glueball states on its trajectory. This trajectory should have the observed rather small 



slope (Jp and "supercritical" intercept ap(0) > 1. The calculations done so far, while not complete, were 
at least able to confirm this concept of the Pomeron. 

Due to the smallness of the gauge coupling, the two gluon state is expected to dominate the Pomeron 
in reactions involving high transverse momenta. However, at high energies or correspondingly for high 
rapidities, a simple perturbative calculation of the two gluon exchange is not sufficient and interactions 
between the two gluons must be taken into account. This is done in the well-known BFKL Pomeron 
in leading logarithmic approximation (LLA) in log(l/x) with a fixed QCD coupling constant. This "hard" 
Pomeron can be represented by the exchanges of gluon ladders. 

An interesting and characteristic feature of the BFKL Pomeron is the diffusion in the transverse 
momenta (in fact in the logarithms) of the gluons as the rapidity increases along the gluonic ladder, ft can 
lead to a growth of the mean transverse momentum with energy and one expects that at very high energies 
the "hard Pomeron" dominates in processes with high transverse masses. In this region of high virtualities 
also the DGLAP evolution might be applied. For the diffusion property at fixed as it is essential that the 
BFKL equation is scale invariant. Clearly it will be interesting to examine how the situation will change after 
this invariance is broken. Since the main source of the scale invariance violation is the dependence of the 
effective gauge coupling on the QCD scale A, the most obvious way to study this question is to include the 
running gauge coupling in the BFKL equation. Evidently there are other higher order corrections: one part 
concerns the BFKL gluon ladder itself , but one expects that they should not change its qualitative features 
much more; other corrections involve multiple Reggeon diagrams Q| . Some aspects of these questions have 
already been considered in Refs. [5[ |TT|. 

It is the aim of this paper to analyse the generalization of the BFKL approach to a running coupling 
constant in the most straightforward fashion and to arrive at a comprehensive physical interpretation. Our 
numerical study in Ref. |12| of the BFKL equation with a running coupling constant has already shown that 
the BFKL Pomeron is of purely perturbative character only at the ends of the gluon ladder, while becoming 
progressively soft towards the middle. We study here the presence of a running coupling constant in more 
detail in an analytical fashion and, in addition, examine the singularity structure of the Pomeron in the 
complex angular momentum plane. 

The outline of our paper is as follows. We begin with a review of the original BFKL equation with 
fixed coupling constant in Section 2. Then in Section 3 we present the exact solution of the BFKL equation 
with running a^, chosen in the simplest form for the full range of transverse momenta and without an 
infrared cut-off. The main result is that we find an essential singularity in the complex angular momentum 
plane, independent of any transverse momentum scale, that leads to an energy dependence for the cross 
sections of the form e'^v^. This result is further explored in the following sections. In Section 4 we examine 
a diffusion-like approximation to the full operator equation. We find that this diffusion solution contains 
an infinite sequence of poles in the complex variable the conjugate of the rapidity y. This resembles the 
series of poles found earlier by Lipatov which replaces the cut in the case of fixed ag. We examine the 
residues of these poles and their importance in hard processes. A simple quantum-mechanical analogue is 
presented to make the essential characteristics of the diffusion solution more transparent. In Section 5, the 
contribution from the singular part of the BFKL kernel is discussed. It represents large increments in the 
transverse momenta of the gluons along the ladder and is shown to lead to a DGLAP-like behaviour. Based 
on the features we found for the BFKL Pomeron with running as , we arrive at a simple more general model 
for the Pomeron; it reflects all the essential properties of the exact solution of Section 3. Its main point is 
that the Pomeron is of soft, nonperturbative origin, but can have hard ends, when probed by a hard device 
as in deep-inelastic scattering. In the limit of very large rapidities the hard ends become small compared 
to the soft part. Applications to onium-onium scattering are considered. A summary and our conclusions 
are contained in Section 6. The main conclusions of this paper were already reported in Ref. p3|. 
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2. BFKL EQUATION WITH FIXED COUPLING CONSTANT 



In this section we review shortly the BFKL equation with fixed coupling constant ag. It has been 
experimentally established that the hard cross sections corresponding to the exchange of vacuum quantum 
numbers have a power-law increase as function of the CMS energy ^/s according to cr ~ ^A^^p with Aexp 
0.3. A simple two-gluon exchange as model for the hard Pomeron g however leads to constant total cross 
sections, and diagrams leading to large log(s) contributions must be taken into account. BFKL included 
such more complex diagrams which could be resummed and represented as exchanges of effective gluon 
ladders. In these ladders non-local gauge invariant Lipatov vertices and reggeized gluon propagators are 
the building blocks (see Fig.l). This was done under the assumption of "multiregge kinematics", where the 
rapidity monotonically increases along the ladder in large steps, 5^1. These ladder-like configurations 
make the cross sections increase fast with energy. BFKL performed their calculations for fixed coupling 
constant as in leading logarithmic approximation, i.e. 

as < 1 ; UsV^l, (2.1) 



where the rapidity interval between the ends of the ladder, y, is related to the total energy, according to 
y log(s). The simple ladder structure of the diagrams makes it possible to write a Bethe-Salpeter type 
of equation for this amplitude. This leads to a linear integro-differential equation, known as the BFKL 
equation for the low x behaviour of the unintegrated gluon distribution in the nucleon: 

dfjy, fcj) ^ r°°dql 



/(y,<?i)-/(y,fci) , /(y,fci) 



kl\ 



asC{kl)(^f{y,kj 



(2.2) 



where we defined the rapidity as y = los(xo/x) with xq a higher x value at which the initial condition 
should be given, and as = 3as/7r. In Eq.( |2.2[) the integration over the angle between k± and q± has already 
been performed. The function f{y,k'^) is proportional to the imaginary part of the virtual gluon forward 
scattering amplitude with vacuum quantum numbers in the t-channel; the first term of the integrand 
corresponds to the emission of real gluons while the second term originates from the virtual part of the 
amplitude. In a reaction with characteristic scale Q^, one probes the integrated gluon distribution, 

Giy,Q')= r ^f{y,kl) . (2.3) 



Since the BFKL equation is scale invariant, it can be solved in a simple way using Mellin and Laplace 
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transformations with respect to A:^ and y. Therefore we define 



oo 



dfc2 / fc2 \ 



/(?/'/3) = y^ /(?/,fci)=y_ due-^-f{y,u) , (2.4) 

f{Lo,u)= / dye-^yf{y,u) , (2.5) 
Jo 

/•oo 

F[u,p)^ / dye~^yf{y,f3) , (2.6) 



where u = log(fc^/m2), a definition we use throughout the paper, and where m is a scale introduced in 
order to make the transformed distribution /(y,/3) in Eq.(2.4) dimcnsionless. The variable uj, conjugate to 
the total rapidity interval y, has a physical interp reta tion relating it to the (complex) angular momentum: 
Lo = i — \. Applying both transformations to Eq.(2.2) leads to the algebraic equation 



uF{uj, 13) = /(O, /3) -t- asC{l3)F{uj, /?) 



(2.7) 
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with the solution 



/(0,/3) 



(2.8) 



where 



= 2vl/(l) - _ v|/(i _ /3) , 



d\ogV{l3) 
dp 



(2.9) 



is the eigenvalue spectrum of the BFKL integral operator in Eq.(2.S). The Mellin transformation with 
respect to k± was only well defined for < Re /3 < 1 . In this range of (3 the function £ has a minimum at 
P ~ 1/2 while it is singular at the edges P — and P — 1 (see Fig. 2). The limit P ^ Q corresponds to the 
hard limit, where u is large. 



We take the inverse transformations over P and lu of Eq.(2.8) to get back to f{y, u) 

pc-\-ioo 



2m 
1 

27ri ./d 



dcj e"^F(w,/3) 



C — 200 



d/3e"''/(y,/3) 



(2.10) 
(2.11) 



with c a real number to the right of all singularities in lo and d an appropriately chosen real number between 
and 1. This yields the solution of the BFKL equation 



/(2/,u) 



dv! f{0,u')gf{y,u-u') 



(2.12) 



where /(O, u') is the boundary condition at y = 0, and with the Green's function 



2-111 



dp exp {j3u + yasC{P) 



(2.13) 



where we took d = 1/2 due to the pinch of the P contour at the point P = 1/2. This pinch results in 
a leading singularity in Eg. (pT^) at uj = asC{\/2). For large y the saddle point method can be used to 
calculate explicitly the Green's function. With the help of the expansion of the kernel around P = 1/2, 



a,£(/3) « A + i?(/3 - 1)2 



A = 4a,log(2) , S = 14a,C(3) , 



(2.14) 



one then obtains 



^jAiiBye-^ 



exp yA 



(2.15) 



The character of the corresponding singularity in lj and its dependence on u can be more easily seen in the 
(w, M)-representation where the Green's function takes the form 



1 



^4Be-"(w - A) 



exp 



B 



(2.16) 



which shows the branch point singularity at a; = A. 
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From Eqs.(|2.12D and we see the two characteristic features of the BFKL Pomeron. First, the 



Pomeron intercept A ~ 2.6 is relatively large, of order 0.5 for a realistic as ^ 0.2. Therefore the observed 
energy dependence of cross sections by the HERA experiments, a ^ s^""'', was seen as a possible sign of 
the appearance of the BFKL Pomeron in hard scatterings. Second, the diffusion of the u distribution as y 
grows: < {u — m)^ > ~ 4i?y. As will be seen later, this important diffusion property is a consequence of 
scale invariance when ag is fixed. 

We conclude this section with two comments about the applicability of the BFKL Pomeron in LLA. 
Due to the diffusion of the u distribution, it can happen that as y — > oo one reaches regions of small 
transverse momenta, where non-perturbative effects take place. For the perturbative BFKL Pomeron to 
be applicable, these contributions should be small. One can simply exclude transverse momenta below 
some fixed fc™*" in all gluon ladder diagrams by introducing a cut-off in the integral that enters into the 
BFKL equation [p^jjl^ ]. Another way to achieve this is to consider the BFKL Pomeron in onium-onium 
scattering. In this case there is a purely perturbative regime near both ends of the gluon chain due to the 
small transverse sizes of the onia, R\ ~ 1/^1- They can be chosen sufficiently small, so that the transverse 
momenta in the gluon ladder do not reach the infrared region in the entire rapidity interval between the 
two onia. In the numerical examples in Ref. [|l2|, the influence of the small fc^ region on the gluon ladder 
was found to be not large for fixed as', the giuons in the ladder are in the mean far from the soft region. 
However, for running as the diffusion to the infrared region was shown to proceed much faster. 

Concerning the assumption of multiregge kinematics, we can make the following estimate for the 
BFKL Pomeron, based on a general property, independent from the underlying QCD dynamics. For a 
ladder diagram with vector particle exchange the cross section grows as e^^. Expanding this exponential 
function as X]n(2^^)"/"'i 0^*2 obtains the cross section in terms of contributions from ladders with n rungs 
and a corresponding mean multiplicity n — yA. As result, the mean rapidity interval between neighbouring 
rungs is 5 =< yt — yi+i >~ A^^. This quantity is not large: for as — 0.2 one has 5 ^ 2 or A ~ 0.5. 
However, the BFKL was derived under the assumption of multiregge kinematics, i5 3> 1. If we consider 
ladders making up the Pomeron with a lower cut-off, i5o — 1, for the rapidity difference 5 = yi — yi+i, a 
sizeable fraction of the giuons are eliminated. This is due to the fact that the giuons are concentrated at 
small rapidity differences as was confirmed numerically in Ref. For as — 0.2 it was found that the new 
A characterizing the resulting energy dependence of the cross section becomes smaller by approximately a 
factor of two. 



3. EXACT SOLUTIONS OF BFKL WITH RUNNING as 



In this section we include a running coupling constant in the BFKL equation. The procedure to 
do this is not unique. It is often required that the so-called "bootstrap condition" ||^ is obeyed. This 
consistency condition follows from the requirement that the solution of the BFKL equation in the gluon 
colour octet channel should be that of the reggeized gluon itself. Further one should check that in the limit 
of very high transverse momenta the double logarithmic approximation of the DGLAP equation is obtained. 
Nevertheless, many ways to introduce a running gauge coupling ag in the BFKL equation remain possible, 
but they are expected only to result in non-essential corrections. 

We include a running coupling constant in the BFKL equation by using instead of the fixed one in 
Eq.( p.2| ) the perturbative expression 

"-^'^) = 51og(4/A^) ' ^'-'^ 

where A is the QCD scale parameter and b = (33 — 2iV/)/127r. This choice was also made in Ref. ||^,^. We 
put the running coupling constant under the integral in the following way: 



df{y,kl) ^ ^ r dqi^2 
dy Trio <li ^ 



s{ql)f{y,ql)-as{kl)f{y,kl) , a,(fci)/(y, fci) 



(3.2) 
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This has precisely the structure of the original BFKL equation when we write h{y, fc^) = as{k'^)f{y, fc^): 



dh{y,kj 
dy 



as{kl)C{kl)(g>h{y,kl) 



(3.3) 



It is also possible to leave the running cou pling constant as(fc^) outside the integral. In this case we obtain 
a modified BFKL equation as in Eq.(3^), but now for f{y,k\) instead of h{y,k\). Therefore it yields a 
gluon distribution with an additional factor as(k\) as compared to our choice where it is put under the 
integral. As we will see later, our method is the one that leads to the correct behaviour for large transverse 
momenta. 



We begin with a qualitative argument in order to show how including a running in the BFKL 
equation can modify the characteristic BFKL features. For constant the large rapidity behaviour of the 
total cross section predicted by the BFKL equation is given by (j{y) ^ exp(?/asl2 log(2)/7r). The as entering 
this expression should be taken at some fixed, effective value of fcj^ that must be of the order of the mean fc_L 
in the giuonic ladder representing the Pomcron. However, for moderate transverse momenta at the lower 
end of the gluon ladder and a high virtuality at the top of the ladder {e.g. in deep inelastic scattering), 
the mean \og{k\) changes along the ladder and is proportional to Therefore, a simple estimate for 
the effect of a running coupling constant can be obtained by using ~ 1/ y^, leading to a cross section 
<j{y) ~ e'^v^. This behaviour is also what we find in the double logarithmic approximation of the DGLAP 
equation. As we will see later, in the complex angular momentum plane lo ~ j — 1 it corresponds to an 
essential singularity of exp (l/w) type. 



By using as, Eq.(3.1), for all k±_ we partly suppress the infrared region, which is in line with our goal 



to investigate the influence of a running ag on the structure of the BFKL Pomeron at large /cj_ . Qualitatively 
this may be understood from the fact that the contribution of a ladder for vector exchange, with n rungs 
and in a limited region of integration over fc^ , results in amplitudes of the form 



with as{^i^)'D{^i, ^12) 



d?k^ a.(fci) 



(3.4) 



where p, the mean scale at which the running coupling constant has to be taken. In this way one obtains a 
Pomeron intercept with Aq ^ as{il'^)'D. Hence we see that the simplest way to suppress the soft contribution 
is to make ag negative for small fc_L so that the corresponding sing ularity moves to the left half w-plane. 
Therefore it is not bad that the standard expression for a^, Eq.(3.1), has this property for k\ < h? . Note 
also that the singularity in Eq.( |3.1| ) is integrable at k\ = . 



As in the case of constant as, we solve the BFKL equation with running as, Eq.(3.2), by using the 
Laplace and Mellin transformations of Eqs.(2.4) and (2.6). In the McUin tranformation of Eq.(2.4) we set 



A . With the help of the simple identity 



fci 
A2 



log 



A2 



d(3 I A2 



(3.5) 



we then obtain for the Mellin transformed differential equation for the distribution h{y,u) 

d^h{y,P) 



dydp 



asC{P)h{y,(3) , 



(3.6) 



where C{(3) is the standard BFKL kernel, Eq.(2.9), and we d enoted as — 3/7r&. Further we perform the 
Laplace transform over the rapidity y which yields for Eq.(3.6) 



(3.7) 
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with Hq{I3) a function determined from the initial condition at y = 0, 

dh{0,l3) 



dp 



(3.8) 



The homogeneous solution of Eq.(3.7) is 

Hh{uj,P) = i7if(w,/3o) exp 



(7^(/3) - 7^(/3o)) 



(3.9) 



where 



7^(/3) = / d(3'C{(3') = 2vI/(l)/3 - log 



m 

r(i - p) 



(3.10) 



The function TZ{/3) is plotted in Fig. 2 . Thus the total solution becomes 



(3.11) 



This is the exact solution of the BFKL equation with a running coupling constant. 



It is essential that Eq.(3.11) has only a fixed singularity at w = 0, i.e. not depending on (3. This 
is in contras t to the situation with a fixed as where the analogous solution in the (w, /3)-representation is 
given in Eq.( |2.§| ), and results in a /? dependent u-pole. For running as there are no pinch and endpoint 
singularities when carrying out the inverse tranformations needed to obtain the solution as function of y 
and u, since also the singularities in f3 are not dependent on lu. 

For the purpose of studying the y dependence or, equivalently, the uj structure, we want to focus 
on the first term in Eq.(3.11). The second term does not introduce any new singularities in uj since its 



integrand is proportional to the first term and singularities in uj are independent of /?. The solution of the 
BFKL equation with running coupling constant then becomes in the (w, u)-representation 



ft.(w, u) = h{u!, Uq) 



Kr{uJ, u) 



(3.12) 



with Kr{i-j,u) given by 



Kriu;,u) = — 



d-\-iQo 



dp exp 



d—io 



(3u -n{(3) 

UJ 



(3.13) 



The function ft.(tt', uo) is determined by the boundary condition for u at uq. The function hly, u), which will 
yield the gluon distribution as function o f rap idity and transverse momentum, can finally be obtained from 
the inverse Mellin transformation of Eq. ( ^.12 ) with respect to 



iO. 



We now discuss an explicit solution of the BFKL equation with running coupling constant at very 
hard scales, appropriate for deep inelastic scattering. In this case we expect to find the double logarithmic 
approximation in log(l/a;) and log(A:^) of DGLAP. The large transverse momentum behaviour is determined 
by the low /3 limit of Eq.( 3.13| ). We can then use 



r(/3)]"^/"~^-"=/'" , (3.14) 
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and the inverse MelHn transformation in Eg. ( 3. 13 ) can be easily evaluated. Only the singularity ai (3 '-^ 
is important; contributions to h{uj,u) from cuts starting at /3 = —1,-2,... have additional factors e~", 
e~^", . . . and thus can be neglected at large u. After deformation of the path along the imaginary axis we 
can use Hankel's contour integral formula for F-functions ||l^ and obtain 

Kr[uj,uo) 2711 \(JJ Kr(uj,uo)r{as/uj) 



The interpretation of this solution is more transparent if we rewrite it in the more familiar form 



/(w,w) 



as{u) 



h{uj,u) /(w,uo) 



as{u) 



with 7<,g(w) 



(3.16) 



In this form for /(w, u), we recognize the anomalous dimension jgg{Lu), thus showing that wc indeed arrive 
in this limit at the DGLAP solution. As was remarked earlier, we could have placed the running coupling 
constant outside the integral in the BFKL equation. We would then have obtained as a result of that 
equation directly the function f{uj,u), the gluon distribution. In that case the gluon distribution differs by 
an extra factor as(u) ~ 1/u, yielding a discrepancy from DGLAP. 

An important ingredient in the solutions is the boundary condition at low transverse momenta 



/(w, uo) which provides the function h{uj,uo) in Eq.(3.12). Poles in /(w, uq) at = result in a ra- 
pidity dependence e^**'. A choice for the boundary condition is that there are no poles in the right half 
plane and thus no supercritical contribution originating from the region of small u. The singularities of 
/(w, uq) are now to the left of w = and thus not important, and the hard part of the Pomeron dominates. 



In that case the dominant part in the limit ~ is contained in the function Kr{uJ,u) of Eq.(3.13). Such 
an approach which is only based on the Kr{oJ, u) and not on boundary conditions could in the past provide 
a reasonable description of the first HERA data. The present low x data show the need of a supercritical 
intercept input. The Mellin transformation from lo to y can be performed and yields a modified Bessel 
function Ii. Thus one finds for Kj,(y, u) after inverse transformations over lo and (3 



(3.17) 



For large rapidities one may use Ii{y) ~ c^/ ^/2ny so that Eq.(3.17) simplifies to 

pd+ica / ~ -T) 1 \ 1/4 



^rM-hLZ'^{W) exp(./3 + 2y^^) . (3.18) 



By using the method of steepest descent, the inverse Mellin transformation then yields the standard double 
logarithmic approximation to the DGLAP equation 

Kr{y, u) ^ "gj^2^" j exp (^V45^ylogu) • (3.19) 



This solution corresponds to the well known double logarithmic approximation where not only a strong 
monotonous increase of rapidity is imposed, but also a strong ordering of the transverse momenta along 
the ladder. It thus means that the resulting cross section grows as e'^^, slower than for the observed 
"supercritical Pomeron", which behaves as e^'^""'''. 

The above boundary condition seems not to be realistic if one tries to decribe the soft high energy 
processes, by assuming the existence of only one Pomeron, containing both soft and hard mechanisms. In 
the case where we expect the existence of a supercritical Pomeron (ap(0) — 1 + Ao > 1) coming from QCD 
interactions in the nonperturbative region, we can choose as boundary condition 
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I, \ JO 
/(w,Uo) = T- 

w - Ao 



(3.20) 



Since Aq > is to the right of the essential singularity of Kr{uJ,u), the boundary condition gives the 
dominant contribution at large y while the singularity of Kr{uj, u) at uj = represents the hard corrections. 
(Strictly speaking, at the boundary u = uq the behaviour of the gluon density may have to be chosen also 
in accordance with the BFKL equation because, depending on uq, the perturbative dynamics contained the 
BFKL equation may already apply.) To analyse the influence of the soft, low transverse momentum region 
in the BFKL equation one can introduce a lower cut-off in the integrals over k±^ in the BFKL equation. 
But it is more complicated to apply this method to the case of a running as than for a fixed one ||l4,|l5|. In 
the next sections, we will discuss in more detail the physical picture of the gluon ladder for this cho ice o f 
boundary condition corresponding to supercritical behaviour. If one supplies the boundary condition ( 3.20D 



the mean transverse momenta in inelastic processes will not grow with energy; this will also be seen later. 

The main result of this section is that the hard part of the BFKL Pomeron for the case of a running 
as leads to a cross section behaving as e'^^ , growing slower with y than the usual "supercritical" Pomeron. 
The situation where the leading y dependence is found to be determined through the boundary condition, 
rather than the BFKL kernel, will be discussed below. 



4. THE DIFFUSION APPROXIMATION TO THE BFKL KERNEL WITH RUNNING as 



We now discuss two main aspects which determine the character of a parton cascade represented 
by our effective gluon ladders. They are the branching of the gluons leading to an exponential growth of 
their number with rapidity and the gluon diffusion in the transverse momentum plane. Both determine how 
the cascade evolves from the fragmentation region of a fast particle towards "wee" partons. Of particular 
importance is whether the branching and diffusion velocities depend on k±. In asymptotically free theories 
the branching velocity becomes small for high k±. In this case all properties of the parton distribution at 
large k± will be determined by the relative weight of the diffusion and the DGLAP-evolution of partons 
towards larger k±. The scale invariant BFKL equation with fixed as has the unique property that the 
splitting probabilities of partons do not depend on k±. Therefore the distribution of partons in k± will be 
determined entirely by the initial distribution at the end of the ladder, which for large y is spread out over 
an increasingly larger interval of k± by the diffusion mechanism. In this section we will study these two 
aspects for the BFKL Pomeron for running as. 



A. The diffusion-like part in the operator equation 



Due to its convolution structure, it is possible to rewrite the BFKL equation Eq.(3.3) as 

1 dh{y, u) 3 



as{u) dy TT 



C{(3)hiy,u) ■ (4.1) 



Here f3 is an operator acting in the space of functions depending on u. It follows from the form of the Mellin 
transformation that the operator is conjugated to w: $ — A useful property of the operator kernel C is 
that it satisfies the "shifting relation" 

e™/:(/3 + c)e-™ = £(/?) , (4.2) 

which can be proven by showing that both sides have identical Mellin transforms, C{f3). As a consequence 
the "shifted" functions hc{y,u) — e'^"^'^h{y,u) obey the differential equation 
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1 dhc{y,u) 
as{u) dy 



TT 



c/2)hc{y,u) 



(4.3) 



Bec ause we know the dependence of the BFKL integral operator on the variable /3 through C given by 
Eq.( ^.9| ), it is possibl e to fi nd a diffu sion approximation. The function £{(3) has a minimum at /? = 1/2 as 
can be seen from Eq.(2.14). For Eq.(|4.3|) it means that the first order derivative with respect to u vanishes 
when the argument of C equals 1/2. We therefore choose c = 1 in Eq.(4.3) and expand in /3. We introduce 
the function hi{y,u) = e~^/'^h{y,u) and neglect third and higher order derivatives, which leads to the 
differential equation 



1 dhx 

has{u) dy 



A hi 



dv? 



A = 4a,log(2) , S = 145,C(3) 



(4.4) 



This equation resembles a diffusion equation with branching, but with a diffusion coefficient as{u)bB and 
a branching coefficient as{u)hh.. Here the rapidity y plays the role of time and u that of place. This 
analogy will be useful for understanding the structure of the function h{y,u). The diffusion coefficient can 
be interpreted as a diffusion velocity, determining how fast the gluons diffuse in y. Since it is proportional 
to as{u), the gluons diffuse slowly to large u; the region /? ~ 1/2 corresponds to an evolution where changes 
of u in individual steps are not large. Eq.( |4.4| ) must be supplemented with some boundary conditions at 
small u = uq, representing infrared low k± physics not contained in Eq.(4.4). 



For the case with fixed a^, Eq.(l^) is a standard diffusion equation with branching. 



dy 



B 



dv? 



(4.5) 



where the branching and diffusion coefficient are independent of u. The solution h{y, u) — e"^^ hi{y , u) is 
then given by Eq.( p.l2 ) with the Green's function (2.13). It should therefore be clear from Eq.(4.5) that 
for the diffusion behaviour predicted by BFKL scale invariance is essential. 

We have seen above that the region (3 ^ 1/2 corresponds to an evolution where changes of u in 
individual steps are not large: Su ~ as- On the other hand, large jumps in u are related to the region 
/3 ~ 1/m — > 0, where £(/?) ~ 1//3 and for which 5u ~ a^u. In this case we have to set c = and expand 
around /3 = 0; this will be done in the next section. 



With our choice for as{u), Eq.(3J), the differential equation (iA) becomes 



dhi X , , 5 d^hi 

= Ahi + B 
dy dw' 



(4.6) 



A solution of Eq.( [1.6| ) can be simply found by performing a Laplace transformation with respect to y. For 
the transformed function h(uj,u) it follows that it satisfies the differential equation 



7 / \ XT/ \ n d h(uj, u) 

uuj n[uj, u) = A Mu!, u) + B 

dw' 



(4.7) 



where for simplicity and without loss of generality for the present discussion we took /i(0, u) — 0. By using 
the substitution u ~ z{B /loY/'^ + K/to the differential equation (4.7) transforms into the Airy equation, 
/"(z) — zf{z) = 0, which has the solution 



h(uj, u) = Ai(z) 



1/3 
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where Ai is the Airy function that is f inite at z ^ oo. We note that Ai is related to the modified Bessel 
function Ki/3 through Ai(z) (l/7r)y(I73)Ki/3 (2^2/3/3). 



Before we give a full analysis of the solution ( |4.8| ) in the next section, we first examine the dominant 
behaviour at large u. If we assume a boundary condition for hi(y,u) at small u of the form 



(4.9) 



then for I one can try to find a solution of Eq.(4.6) in the form 

hi{y, u) ~ exp (Aoy - Cqu" - du"^^ - ■ ■ ■) 



(4.10) 



Upon inserting Eg. ( 4. 10 ) into Eq.( fl.(]| ) and comparing the leading terms in u which are proportional to 
uhi{y,u), we find that s — 3/2 and A = BCqS^. Furthermore, comparison of terms proportional to 
hi{y,u) gives A = 2s(l — s)BCoCi. Solving for Cq and Ci results in 



= 3/2 , Co^±: 



2 fAo 
B 



1/2 



Ci 



B \ B 



-1/2 



(4.11) 



We have two solutions, positive and negative. Because a fast increase in u must be rejected, only the 
positive solution is accepted for Co 0. An important point is that the hard Pomeron disappears in the case 
of running in this approximation, because the residue belonging to the pole at a; = Aq rapidly decreases 
proportional to exp (— Cqu"^/^) on a scale u ^ ^^"^ ~ 1. (The scale for the disappearance of this solution 
is set by Aq; the closer it is to the critical value Aq = 0, the longer the hard part of the Pomeron survives.) 
A significant large u contribution can therefore only come from the region /3 — !■ in C{f3), which will be 
discussed later. 



The result ( |4.10| ) can be explained in a simple fashion that shows how the transverse momentum in 
the gluon cascade evolves as function of rapidity y in the limit m 3> 1. Assume that starting from the soft 
(i.e. small u) end of the gluon ladder the gluons are in the strong coupling region up to rapidities y — yi \ in 
this region the gluon distribution should then be proportional to hsoft{y — yi,u) ~ e'^"'-"""^'. After that, in 
the remaining part of the rapidity interval yi , the gluons reach large virtualities u o n th e other end through 
the diffusion-like mechanism. This last stage is described by the solution of Eq. (|4.6|) , which for large u 
becomes 



rdiy,u) ~ exp(Ay/u-u^/9B?/) 



(4.12) 



Using this function, we can obtain the gluon distribution for the full y range through the integral 

hi{y,u)^ r dyie'^«^y'y''>hhard{yi,u) . 
Jo 



(4.13) 



This integral can be estimated at large y by the saddle point method which gives 

1/2-1 



hi{y,u) ~ exp 



Aoy - w 



3/2 2 ( Ao 



3 V B 



(4.14) 



^It is interesting that if QCD was not asymptotically free (for example due to a large number of fermions), then 
the diffusion in u described by Eq.(4^) would go in the opposite direction of large u. In such a situation only the 
theory with a fixed limit makes sense when as{u) — > Umax for u — > oo. So in this case the diffusion in u will move 
gluons in the gluonic ladder to higher u and thus the Pomeron will develop according to the BFKL solution (^.15 
with constant ag = amax- 
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coinciding with the main term in Eq.(4.10). In Eq.(4.13) the mean yi is of the order 

1 



< 2/1 > ~ u 



3/2 



(4.15) 



This shows that at large y the hard part of the Pomeron is concentrated on its end. But these st atem ents 
about the structure of the end part of the gluon ladder given by the expressions Eg. ( 4. 10 ) and Eg. ( 4. 15 ) do 
not concern the dominant term, since its residue is very small. 



B. The infinite sequence of poles 



In this subsection we further study the solution (4.8) of the diffusion eguation, not restricting ourselves 



to large u. We pay special attention to the pole structure in the variable uj, related to the angular momentum 
through j — Lo — 1, and the influence of the boundary conditions. Our approach, which seems natural for 
the BFKL and the DGLAP cguations, is to consider solutions only at u > uq, when uo is not small but also 
not very far from the infrared region. On the line u = uq a boundary condition representing the soft physics 
must be given. Of course, it is not possible to choose a uq where only soft physics enter at the boundary; 
depending onu = uq a varying amount of hard physics will already be contained in the boundary condition. 



In the variables uj and u, we thus obtain as solution of Eg. (4. 7) 

Ai(z) 



hi {uj, u) = hi {uj, Uq) 



Ai(zo) ' 



(4.16) 



where /ii(w,Mo) is the boundary condition. The arguments of the Airy functions z = z{u),u) and zq = 
z{uj,Uo) depend on uj and u as given in Eg. (4.8). From the point of view of Regge phenomenology it is 
obvious to assume that hi{uj, uq) is represented as a sum of Regge poles, some of them having a supercritical 
intercept. 



hi{uj,uo) = ^ 



UJ ■ 



(4.17) 



As said above, in ad dition to singularities coming from soft physics, some poles contained in the boundary 
condition. Eg. (4. 17), are connected to the hard physics contained in the BFKL diffusion kernel. 



Poles in Eg. ( 4. 16 ) can also arise fro m ze ros of the denominator. These additional poles occur when 

is negative. This provides an upper bound on the spectrum 



the argument z of the Airy function in Eg. (4 



of the singularities in uj arising from the Airy function. 



A 



i^max < = 4as(uo) log(2) . 

Uo 



(4.18) 



We see that the most right singularity in the complex angular momentum plane, j, is smaller than the BFKL 
intercept calculated with the coupling constant at the boundary scale uq. The position of the zeros of the 
Airy function are with high precision given by the formula ||l^ z = — (37rn/2 — Stt/S)'^^^ with n a positive 
integer. We can find an approximate value for ujmax by inserting uJmax = A/uq — e into the expression for 
z{uj, u). Solving for e, we find 



UJ„ 



A 



1 

2/3 



2/3' 



(4.19) 
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with K = (A/i?)^/^ « 0.47. It is also possible to obtain the locations of the poles in the limit lo A/uq 
where zq ~ — (kA/cj)^/'^. The zeros in the denominator then result in poles located at 

(4.20) 



This shows that in the diffusion approximation with running coupling constant, the solution contains an 
infinite set of poles accumulating at w = for n — > oo. Such a sequence of poles was already found in the 
A(/)|-model in Ref. In QCD this sequence of poles was first found by Lipatov starting from the 

same generalization of BFKL to a running gauge coupling, but using different techniques. The discussion 
of these poles in Ref. |^ is more close to ours. 

To clarify the importance of the poles w^, we calculate their contribution to /ii(y, u). Performing the 
inverse Laplace transformation over uj we obtain 

oo 

/ii(y,M) = ^e^'^'A(A„M) + ^e^'^"AK,M) , (4.21) 

i n—1 



LOn — Ak 



3 3 
-7rn — -TT 
2 8 



where the first term is the contribution from the (mainly soft) singularities A,j of Eq.(4.17) and the second 
that from the Airy function in the denominator; A(r/, u) denotes the residue of a pole a.i lo = rj. The residues 
of the second term, A(a;„,u), initially oscillate in u and can be large up to u '-^ 37rn/2K. After that, for 
uq and large n, they are proportional to 



A(w„,u) 



Ai 



37m 



1/3 



37m 



2/3N 



exp 



2ku 



3/2N 



(4.22) 



We see that for large n, where oJn — > 0, the residues increase and they become "harder". Because the 
residues go quickly to zero, only the largest n values are essential at very high u. The u dependence of the 
residue of a pole at w„ scales with KT^m}^^ . Thus singularities in the sequence a;„ more to the right {i.e. to 
larger values) are softer. 



For large u, it is simple to estimate the contribution from large n to hi{y,u), Eq.(4.21) 

1 



hi{y,u) 



0n(") exp 



2kw\^/^ 2Ak' 



37m 



(4.23) 



where 0n(u) are weakly dependent on u. The largest term in this series is at the point n = (BA^/ttk) jv? . 
Hence the sum can be estimated as 

/ii(y, u) ~ exp [yKju - u^Kj^Bi^ , (4.24) 



which is similar to the solution, Eq.(4.12). So we see that the sum of poles from the BFKL kernel at high 
u can be represented by an effective u dependent intercept with Aeyy(u) = A/u. 

We now comment briefly on the connection between the boundary condition, ft.i(cij, mq) and the BFKL 
kernel in Eq.(4.8). As can be seen from Eq.( [4.18| ), choosing a larger starting value mq moves the upper limit 



ijJmax of the spectrum of poles in lo to the left. In order to lead to the same result for hi{uj, u), independent 
of Mo, this change in the kernel due to the poles near uomax, has to be compensated by a concomitant change 
in the boundary function hi{ijj, uq). For example, a pole at tok can be removed by a zero of hiiojk, uq). If, 
on the other hand, Uq is decreased, a pole A^ originally in the boundary function will move into the BFKL 
kernel. This connection shows that for a stable choice of uq, the boundary has to be in a region where the 
BFKL equation, and especially the diffusion approximation, is already applicable. 
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While we have found an infinite set of poles accumulating towards w = 0, for a fixed coupling constant 
there is a cut in the w-plane (see Fig 3.). It was remarked in Ref. ||] that the sequence of poles w„ for running 
as in some sense corresponds to a cut in the w-plane for fixed as ■ But this analogue is incomplete as far as 
we can see for our case with t — (forward scattering). For fixed as the right part of the cut corresponds 
to the hard Pomeron with mean < u >~ while for running as the right poles are relatively soft with 
< M 1, and only the left side of the sequence uj^ towards uj = becomes hard. These poles may even be 
artifacts due to our particular choice of running as, which allowed an analytical solution. Furthermore, they 
may disappear when other higher order corrections to the Pomeron, for example the mixing with multiguon 
states in the t-channel, are taken into account. These questions go beyond the scope of our paper and need 
more investigation. 

We can get a slightly different view on the evolution of the gluon ladder and the structure of the 
a;„-poles by considering a quantum mechanical analogy. Eq.(4.4) can be seen as the Schrodinger equation 
for a one-dimensional motion in the coordinate u in the interval uq < u < oo. The energy of this motion — 
when time is identified with t = —iy — is given by 



E = 



1 



2m{u) 



f + Viu) 



(4.25) 



where we defined the momentum, potential and coordinate dependent mass, respectively, as 

p =-i(3^ ~id/du , (4.26) 

V{u) = -A/u , (4.27) 

m{u) = u/2B . (4.28) 



In the case of a constant as we have free motion. So when we put the "particle", represented by a wave 
packet, at an initial time at s ome position ui, then the situation at a later time is described by the Green's 
function for free motion, Eq.( 2.15 ), yielding propagation and spreading of the wave packet. But for running 
as the motion is not free; we have a long ranged attractive force in the direction to smaller "distances" 
u. In this potential an infinite "Coulomb-like" series of bound states exist. For the high lying states, the 
particle is located on the average at large distances uq where the motion is quasiclassical. To find their 
energy spectrum, _E„, we use the standard quasiclassical quantization condition 



dup{u) = 



l'A/E„ 
Juo \ 




A - EnU 


2 


(i) 




B 







(4.29) 



which leads to En = (1/n)- (2A^/^/37ri3^/^), coinc iding with Eq.(4.20) in the limit of large n. It is interesting 
to note that the behaviour of the residues in Eq.( 4.22| ) corresponds to the motion of our ficticious particle 
in the classically forbidden region u > A/ En- The momenta in this region are not constant due to the u 
dependence of the mass. 



Pn{u) = \J-2m{u)En = i\ 

A-K n 



2k3 u 



(4.30) 



Using the wave function A„(w) '-^ exp {iupn{u)), we obtain the same result as in Eq.( 4.2^ ) 



So far we have excluded the particle from the region u < uq, corresponding to an infinite potential 
barrier. We can change the potential (and thus the boundary condition at uo) in order to allow the "particle" 
to go to smaller u. This can for example be achieved if we choose the potential 



V{u) = 



log {v{u) + e") 



(4.31) 
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where the function v{u) — > at w — > cx) and v{u) —fVo>latu—f —oo. At large u > we have the same 
potential as before but now we can put uq = —oo . In this case the motion in m < will be unlimited also 
for some E < 0, and all discrete eigenvalues disappear. In addition, we can choose v{u) in such a way 
that V{u) will have a local minimum at u ^ 1 so that we generate a low lying bound state. In terms of our 
pole structure in to, this corresponds to a "soft" state with a high intercept resulting from this potential. 

In Ref. (19) the structure of the BFKL Pomeron for t < was studied. It is amusing that we can 
arrive at their findings concerning the Pomeron trajectory by using the following shortcut based on our 
quantum-mechanical analogue. The position of the Pomeron pole at some momentum transfer t < can 
be simulated from our t = example by restricting the range of the "coordina te" u {i.e. the range of the 



transverse momenta, through the condition u > uq — log {~t). We choose in Eq.( [4.29[ ) a large uq and simply 
put n = 1 to obtain the leading BFKL Pomeron trajectory, 

u;p{t) E,^ ^ , (4.32) 

uo{t) 



where uq = l og j— t/A'^) ^ 1. The trajectories for the other "satellite" poles also can be estimated for large 
uo from Eq.([4.29[), 



w„(uo) ^ A(uo) 



2k Uq 



2/3' 



(4.33) 



These relations for the Pomeron trajectory coincide with the ones that were found in Ref. |l9| in a precise 
(and therefore more complicated) way. The above results also show the explicit dependence on the infrared 
cut-off Uq. 

In conclusion, we repeat that the general significance of the poles at w„ needs future investigation. 
They may be a feature closely connected to our particular choice of the running coupling contant. It will 
be shown below that for large u their contribution is small in comparison with the /3 ~ part of the kernel, 
corresponding to the DGLAP mechanism. 



5. POMERON WITH HARD ENDS 



Eq.(ll) has the structure of the kinetic equation. The differential equation (4^) discussed above 
corresponds to the diffusion approximation, where fast increases in the "coordinate" u do not occur. For 
the BFKL Pomeron with fixed as such jumps in k± are not essential, because the splitting of gluons is 
universal for all k± and the diffusion mechanism is sufficient to evolve the gluon distribution in the "time" 
y to obtain f{y,u) at larger u. All possible inhomogeneities, arising from initial conditions or the infrared 
region, are smeared out by diffusion. 

As we saw above, for the BFKL Pomeron with running ag the situation is different. In the diffusion 
approximation, for which the vicinity around (3 — 1/2 of £ is essential, the gluon splitting probability 
Bas(u) at large k± decreases with u. Therefore only a small fraction of the initially moderate k± gluons 
diffuses to the region of larger k±. As result the density f{y, u) ~ ue~^/'^h{y, u) decreases fast with u. This 
has been discussed in detail above. 

In contrast to the diffusion-like behaviour resulting from the expansion of the BFKL kernel C around 



its minimum, the limit /3 — > of the kernel leads to a growth of /(y, u) as function of u as shown in Eq.(3.19) 



For a BFKL Pomeron with running Us only this part can generate a significant contribution at large u. In 



this section we briefly show how these previous results, the low x DGLAP solution, Eq.(3.19), follow quite 



directly from the BFKL operator equation Eq.(4.1). 
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It is well known that a DGLAP-like evolution equation follows simply from BFKL when one takes 
into account only configurations strongly ordered in fcj^; this limiting case corresponds to the part of the 
eigenvalue spectrum C{l3) situated at /3 ^ 0, where C{P) ~ 1//3. In the operator equation (11), the rapidly 
increasing part of the BFKL Pomeron is found by using the approximation to the BFKL operator 



L0) c r 1 = 



du I 



(5.1) 



Thus instead of Eq.(4.4) we get the integral equation 

1 dhi{y,u) 



bas(u) 



dy 



du hi{y,u) + hi{y,ua) 



(5.2) 



This approximation is plausible if one looks at Eq.( |2.2| ): take ^ q±_ in the integral of E g. (4.1), which 
then transforms directly into Eq.(5.2) if we introduce the same integration limits. Eq.(5.2) can be solved 
by means of a Laplace transformation in y and taking the derivative with respect to u, which yields 



hi{y,u) 



c-\~ioc 



duje'^yhi{uj,uo) 



.(^o) Y 
's{u) J 



(5.3) 



which is the same DGLAP asymptotic solution as in Eq.( 3.15| ) 



From our considerations follows a simple model for the Pomeron suitable for cases when one probes 



the Pomeron in a hard process, as in deep-inelastic scattering. In fact it is prompted by Eq.(5.3), where 
hi{y,u) is written as the inverse Laplace transformation of the product of two functions: the first function 
is the boundary condition function /ii(w,wo), containing mainly the soft, non-perturbative physics, and 
the second one is the BFKL kernel with running as, Khard{^,u), containing the essential singularity in 
u). For the boundary condition we expect a supercritical behaviour and therefore use the function given 



in Eq.( 3.2C ). We use this structure to directly model the unintegrated gluon distribution f{y,u). For our 



model we choose ^ 

f{y,u) = f„l dy,e^y-y'^^°Khard{yi,u) , (5.4) 



y 







where 



Khard{y,u) ^ — / dLOe'^yKhard{uJ,u) . (5.5) 
-^T^l J c—ioo 

Since we are interested in a hard process where u is large, we can model Kfiard{y,u) by a DGLAP-type 
kernel that contains all essential features found in Section 3: 



Khardiy,u)^Io{2y/yau)) , (5.6) 

with Jo the zeroth order modified Bessel function depending on the "invariant charge" 

^(w) = a,log — -— . (5.7 
as(u) 



■^A similar expression was already proposed in Ref. in the context of the linked dipole chain model In 
Ref. it represents an interpolation between DGLAP and the hard BFKL Pomeron. We thank G.Gustafson for 
his remarks. 
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In Eq.(5.4). the transverse momenta of the gluons are large up to yi. The values of yi where the dominant 
contribution to the integral in Eq.(5.4) comes from, can in general depend on the total rapidity y. However, 
since the y dependence of the integrand factorizes in our model, it is obvious that the dominant yi region 
at very large y only depends on u, being independent of y. We calculate the dominant interval yi below. 



Expanding the modified Bessel function Iq, one obtains 



1 °° 1 

/(y,«) = /oe^^'' — ^- 



n=0 



7(n+ l,yAo) 



An 



(5.8) 



where 7(n, x) is the incomplete F- function. 

For large y, the Bessel function can be approximated as 



Ioi2^/y^)) 



1 



exp 



(5.9) 



and the integral in our model, Eq.(5.4), can be extended to oo. It then can be done analytically 
f{y, u) = fo dyi e^y-y^^^^>Io{2^^J^)) ~ /oe^«^ / dyi e-y^"^" Io{2^^^) 

JQ Jo 

3 

Ao Ao V "0 



(5.10) 
(5.11) 



This answer can also b e ob tained from Eq. (|5.8|) by noting that for y ^ oo the factors between brackets 
approach 1 so that Eq.(5^) reduces to Eq.( ^.ll|) . The solution shows, as expected, a y dependence char- 
acterized by the Pomeron intercept Aq that originates from the soft boundary, and has multiplicative hard 
corrections, which again involve Aq. For large ?/, the dominant contribution to the above integral can be 
easily shown to come from the region around a saddle point at 



Vi 



(5.12) 



The quantity Ao ~ 0.1-0.3 is small, and depends on the chosen boundary condition (i.e. on the intercept of 
the Pomeron at scale uq). The invariant charge ^(u) starts to grow relatively fast at uq as can be seen from 
Eq.( p.7|) . The subsequent growth becomes very slow due to its doubly logarithmic form. For a boundary 
scale~fcg = 2 GeV^ and A = 0.2 GeV, one obtains ^(u) = 0.3 at fc^ = 5 GeV^, ^(u) = 1 at fc-f = 10^ GeV^ 
up to £,{u) = 1.7 at a very large scale, fc^^ = 10^ GeV^. Since the denominator in Eq.(5.12) is small, the 



(s) 

saddle point yi is already large - of order 10 - at moderate transverse momenta. 

A similar estimate for the relevant rapidities yi can be obtained by calculating the mean yi , weighted 
with the normalized integrand in Eq. (5.4): 



< yi >= (e(") + Ao)/A^ ~ y\ 



is) 



(5.13) 



For the fluctuation around the mean of yi one finds 



which shows that the fluctuations in yi around its mean are proportional to 1/ \fy^ and thus small for 
yi"^ » 2/Ao. 
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In current experiments the available maximum rapidity y is of order 10. For such energies the saddle 
point y[^'' can be large with respect to y and even exceed y. As a consequence, it is then not correct to 

exte nd t he integration range in Eq.(5.1C) to infinity. If y/y'^^ is small, a more adequate representation of 
is given by 




/(2/,u) = /o> In{2.M^)) . (5.15) 



Coming back to Eq.(5.4), we see that our model provides us with a prescription how to incorporate 
the hard part of the Pomeron, Khard{y, u), in reaction amplitudes based on the exchange of Reggeons when 
y is large. In this framework it is simple to calculate other quantities. In the amplitude for "onium-onium" 
scattering with virtualities ui and U2^ the contribution of the Pomeron can be written as 



A{y,ui,U2) 



dyidy2dyz Khard{yi,ui) foe^°^^ Khard{y3,U2)S{y - yi 



y2 - ys) 



(5.16) 



It corresponds to a Pomeron that contains soft physics, in between two hard ends and can be diagrammat- 
ically represented by Fig. 4. In the language of Regge theory, this contribution is an enhanced diagram. In 
the limit of ?/ — > cxD the ends can be considered as vertices, 

A{y,u,,U2) ^ g{ui) foe^"y g{u2) , (5.17) 

where the hard vertices have the asymptotic form 

5(m) = <?oexp(C(«)/Ao) . (5.18) 



One can apply the hard probe to the Pomeron not only at the ends. When we calculate the inclusive cross 
section of th e production of a hard gluon (jet) with rapidity in the central region, we get a result analogous 
to Eq.dslel) . If we assume for simplicity that ?/ — > cxd, so that the hard parts become hard vertices as in 
Eq.( [5.17 ), we obtain 



fVi ry-yi 

Anci{y,yi,ui,u,U2) = fog{ui) g{u2)ey^'> / dya dyb Khard{ya,u)e~^°^'"'''^'"'''^ Khard{yb,u) 

Jo Jo 



(5.19) 



where yi and u are the rapidity and \og{k'j^) of the measured gluon(jet), respectively. 

Finally, we compare the fixed and running as cases for onium-onium scattering, which is usually 
considered as a "laboratory" for the perturbative BFKL equation. For simplicity, we choose similar onia 
with small transverse sizes i?^ . The transverse momenta of the gluons at the ends of the ladders are then 
large: ui = U2 ^ log jj^^^- The coupling constant at these scales is of order 1/ui. Because in our example 
the two onia arc equivalent, the internal structure of the Pomeron is symmetric around the middle, where 
the rapidity is half of the large full interval y. Therefore we only consider the first half of the ladder. At the 
beginning of the ladder, for rapidities yi ^ ui, the simple two gluon exchange model for the Pomeron is in 
fact adequate, since there is no need for resumming the asyi contributions, which are small in this region. 
This holds for both fixed and running coupling constant. After this stage, for a fixed coupling constant 
one enters into the true BFKL regime. The mean m of a gluon in the Pomeron ladder at rapidity interval 
yi is first approximately constant and of order ui] due to the diffusion the transverse momenta spread out 
with rapidity, < u— < u >>^~ j/i, but this spreading happens equally towards lower and larger u. Then 
at rapidities j/i ~ u\ the transverse momenta of a substantial amount of gluons can reach small values, 
M ~ 1. A boundary condition supplied at low tranverse momenta then "prohibits" the transverse momenta 
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u to diffuse into the soft region. As a consequence, the mean transverse momenta in the lad der st art to 
grow in this region. For running as the sitation is different; now the amphtude is given by Eq.(p^). The 
transverse momenta, or u, decrease in large steps, but remain large up to rapidities yi ^ ^{u)/Aq. At this 
point they reach the soft regime. In the limit of very large y the ends are relatively small and the soft part 
covers the most part of the rapidity interval y. 

We thus obtained a model for the Pomeron in hard scatterings that lends itself to a simple diagram- 
matic interpretation. In this model the Pomeron with a running as is of soft, nonperturbative nature, but 
it has hard ends when probed in hard scatterings. The hard ends become small in the limit of very large y. 
However, in the current deep-inelastic scattering experiments the total rapidity intervals are to small to see 
the soft nature of the Pomeron since the hard ends are large, thus leading to large hard corrections. 



6. CONCLUSION 



The original derivation of the BFKL Pomeron was carried out for a fixed coupling constant in the 
leading logarithmic approximation, a^y ~ 1. The physical picture behind it is a gluonic ladder where the 
rungs are strongly ordered in rapidity, the so called multiregge kinematics. Given a characteristic scale at 
one end of the ladder, the transverse momenta of the gluons diffuse. The variance of m = log{k'j_) is of 
the order of y, the total rapidity difference along the ladder. The "hard" BFKL Pomeron applies if this 
diffusion doesn't reach into the soft nonperturbative regime. The gluon splitting along the ladder leads to 
the exponential growth of the gluon density with rapidity, characterized by the BFKL intercept A ~ 2.6 as. 
This intercept follows from a cut in the complex angular momentum plane. 

Assuming a standard form for as{k±), we extended the BFKL equation to the case of a running 
coupling constant. We were able to obtain an analytical solution for the gluon distribution, transformed 
to /3 and w, the conjugate variables of u and y. It has an essential singularity at w = 0, corresponding 
to a cross section of the form e^^. As a consequence, the singularities introduced through the boundary 
conditions at some low scale uq are crucial, because they give the dominant contributions if they are to the 
right of the essential singularity. It was checked that our solution is consistent with the DGLAP expression 
at large y, corresponding to low x. 

We expanded our analytical solution for a running coupling constant around (3 — 1/2. This leads to 
a diffusion-like equation, where the diffusion and splitting coefficients now depend on u. This is in contrast 
to the hard BFKL Pomeron with fixed as, where these coefficients are constant. In this regime the changes 
in u along the ladder are relatively small. We found in this approximation an infinite sequence of poles in 
w, accumulating at w = 0. The poles most to the left are the hardest, i.e. most important at large u. They 
can be effectively represented by a u dependent Pomeron intercept, going to for large u. This behaviour 
could be explained by considering the semi-classical limit of a quantum mechanical analogue. We also made 
a simple estimate for the behaviour of the Pomeron trajectory for t 7^ 0. 

In general, the diffusion-like part of the solution was found to yield a relatively small contribution if 
the changes in the transverse momenta are large. The important part is then due to the /? ~ region of the 
BFKL operator, where large changes in u occur. This is similar to the behaviour predicted by the DGLAP 
equation. We proposed a simple model for large u, which reflects all the essential properties of the exact 
solution: a Pomeron which is of nonperturbative origin, but has hard ends when probed by a hard device. 
In the limit of large u and large y a diagrammatic representation was given, which enables one to calculate 
Pomeron exchange amplitudes quite simply. Onium-onium scattering was discussed as an example in this 
framework. 

In our study we included the running of the coupling constant, an effect we consider the most 
obvious source of the breaking of scale invariance of the BFKL equation. This led to properties of the 
resulting Pomeron, which are very different from the characteristics of the perturbative hard Pomeron 
which one obtains with fixed as . There are of course other higher order corrections which must be taken into 
account. Examples are corrections to the Lipatov vertices and more complicated multi-gluon exchanges. 



19 



Such corrections can still be collected in a modified single Pomeron exchange. There are however also 
contributions to the reaction amplitude beyond the single Pomeron exchange. They include, for example, 
the sequential exchange of more than one Pomeron and the "triple Pomeron" interaction. If these multi- 
Pomeron states become very important, for example at superhigh energies where they are needed to unitarize 
reaction amplitudes, the significance of the single Pomeron as a building block will be reduced. However, 
our study has shown that as long as one uses it as a building block, the running of the gauge coupling must 
be included as it changes the character of the Pomeron significantly. 
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FIGURE CAPTIONS 



Fig.l Contribution of a gluonic ladder with n rungs to Pomeron exchange. Lipatov vertices are denoted 
by dots and reggeized gluon propagators by thick vertical gluon lines. 

Fig.2 Solid line: Dashed line: 11(13). 

Fig.3 a.) For fixed ag: branch cut in the w-plane starting from Re a; = A. b.) For running ag: series of 
poles accumulating at a; ^ starting at Re w = cOmax- 

Fig. 4 Diagrammatic representation of Pomeron exchange in onium-onium scattering. Hatched blocks: hard 
ends. Middle: soft part of Pomeron. 
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